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We derive an exact quantum equation of motion for the photon Wigner operator in non- 
commutative QED, which is gauge covariant. In the classical approximation, this reduces to a 
simple transport equation which describes the hard thermal effects in this theory. As an example of 
the effectiveness of this method we show that, to leading order, this equation generates in a direct 
way the Green amplitudes calculated perturbatively in quantum field theory at high temperature. 
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I. INTRODUCTION 

Classical transport equations for point particles, in the 
absence of collisions (scattering) , have been well studied 
in the past. In more recent years, they have also been of 
considerable use in the study of hot QCD plasma [| . It is 
known, in particular, that in QCD the leading behavior 
of the n-point gluon functions at temperatures T 3> p, 
where p represents a typical external momentum, is pro- 
portional to T 2 and these leading contributions to the 
amplitudes are all gauge independent 0,0,0] . In order to 
extract the leading order contributions to the amplitude 
leading to gauge invariant results for physical quantities, 
it is necessary to perform a resummation of hard thermal 
loops, which are defined by 



P < k ~ T, 



(1) 



where k denotes a characteristic internal loop momen- 
tum. Such a procedure, however, is quite technical and 
the classical transport equation has provided a much sim- 
pler method for deriving the gauge invariant amplitudes 
as well as the effective action which incorporates all the 
effects of the hard thermal loops [alaSBISll3- ^ n sucn 
an approach, one pictures the constituents of the plasma 
as classical particles carrying color charge and interacting 
in a self-consistent manner. 

The basic idea behind the transport equation is to de- 
termine the evolution equation for the distribution func- 
tion. There are basically two equivalent ways of doing 
this. In the first approach one pictures the thermal par- 
ticles, moving in an internal loop, as classical particles in 
equilibrium in the hot plasma whose dynamics are gov- 
erned by classical point particle equations. The transport 
equation can, of course, be derived in a straightforward 
manner once we know the dynamical equations for such a 
particle in the background of a gauge field. For example, 
let us assume that the equations of motion for a particle 
in the presence of a background field is given by 



m— — = k 
dr 



dk a 



eX° 



(2a) 



(2b) 



where r denotes the proper time of the particle and X a 
represents the force it feels in the presence of a back- 
ground gauge field. (In addition, for a particle carry- 
ing color charge in QCD, we have to supplement the 
above equations with the evolution equation for the color 
charge.) The explicit form of X a will, of course, be dif- 
ferent depending on the type of interaction. In general, 
however, the form of X a must be such that k 2 — k^k^ is 
a constant and that the time evolution for fc M transforms 
covariantly under a gauge transformation. Given these 
equations, the classical transport equation for the distri- 
bution function follows in a straightforward manner. 

In the second approach, also known as the Wigner 
function approach, one starts with the Wigner distribu- 
tion function for the quantum field theory of interest in- 
teracting with a background gauge field. The evolution 
equation for the Wigner function is, then, determined di- 
rectly from the equations satisfied by the quantum fields 
aUJl- In the case of self-interacting non-Abelian gauge 
fields, there are two particular issues that need special 
care. First, the Wigner function has to be defined in a 
gauge covariant manner and second, for a self-interacting 
gauge field, one has to use a self-consistent mean field ap- 
proach to identify the background gauge field. We will 
discuss this method more in detail later. But, the advan- 
tage of the second method is that here we do not have to 
know the dynamical equations governing the constituent 
objects moving through a hot plasma. In both these 
approaches, once the transport equation for the distri- 
bution function is known, a current is defined in terms 
of it. This current, which is a functional of the back- 
ground gauge fields, generates the leading hard thermal 
loop amplitudes for the theory of interest. 

More recently, following from developments in string 
theory, there has been an increased interest in quantum 
field theories defined on a non-commutative manifold sat- 

isfymg [a m in in in m m m m ii 



(3) 



where 9^ is assumed to be a constant anti-symmetric 
tensor with the dimensions of length squared. The be- 
havior of hot plasmas in such a non-commutative gauge 
theory and a corresponding transport equation for such 
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theories are interesting questions to study. The difficulty, 
on the other hand, lies in the fact that non-commutative 
field theories inherently describe particles with an ex- 
tended structure [23L 124 . l25| for which classical dynami- 
cal equations are not well understood. In an earlier paper 
[2fil | , we had tried to use the explicit forms of the leading 
hard thermal loop amplitudes in non-commutative QED 
as well as the relation between hard thermal loops and 
classical transport equations to derive a force law for such 
constituents and we had proposed a transport equation 
for non-commutative QED. Such a "phcnomenological" 
equation already exhibits interesting features associated 
with such theories. For example, we have argued that, 
while the charged particles, in such theories, have the ex- 
pected dipole structure, the charge neutral photon field 
exhibits a quadrupole nature in the hard thermal loop ap- 
proximation and that to describe the correct hard ther- 
mal loop amplitudes in the leading order, the classical 
transport equation must necessarily involve "collision" 
terms arising from the extended nature of the constituent 
particles. 

The drawback of our earlier study lies in the fact 
that the transport equation did not have a theoretical 
derivation. As we have alluded to above, this is con- 
nected with the difficulty that we do not understand 
well the classical dynamical equations for extended par- 
ticles in a given background. As a result, in this pa- 
per, we make an attempt to try to derive theoretically 
a classical transport equation based on the Wigner func- 
tion approach, which involves only the properties of the 
non-commutative QED. In Sec II, we describe briefly 
the properties of non-commutative QED and define the 
Wigner distribution function for such a theory. We study 
various properties of this function and derive the trans- 
port equation associated with it. We also define the cur- 
rent associated with such a system. In Sec III, we evalu- 
ate in a direct way the leading order amplitudes following 
from the current, which agree completely with the hard 
thermal loops calculated in perturbative quantum field 
theory. We also compare various features arising in this 
approach with those found earlier and present a brief con- 
clusion in Sec IV. In appendix A, we give some essential 
technical details of the derivation of the transport equa- 
tion. In appendix B, we compile some relations that are 
useful in understanding some other aspects of our calcu- 
lations. 



where c represents the speed of light. Here the Moyal 
bracket is defined to be 



[A,B] MB = A* B - B * A 



(6) 



where the Gronewold-Moyal star product has the form 



A(x)*B{x) 



e i^ a( M n) A(x + OB(x 



C=o=»j 



(7) 



We can, of course, add to this Lagrangian density charged 
matter fields, but since this theory is self-interacting, 
much like Yang-Mills theory, we will continue with this. 
Furthermore, since the photon is charge neutral, this case 
is of more interest, since we normally have transport 
equations for charged particles. We simply note here that 
the theory in Q is invariant under the non-commutative 
U(l) gauge transformation 

A^ix) n- 1 (x)*A ti *n{x) + — n- 1 (x)*d f ,n(x) (8) 

The Wigner distribution function is an operator den- 
sity function in the mixed phase space. For a quantum 
scalar field theory in four space-time dimensions, for ex- 
ample, it is conventionally defined as 



W(x,k) 



-iyk <j) \ [x + l )(t>{x _y ) (9) 



(27rft) 4 

For later convenience, let us define 

i v 

X± = X ± — 



(10) 



The definition in 10 can be easily generalized to the non- 
commutative scalar field theory, by simply introducing 
star products in the product of fields (and treating y as 
a parameter not subject to the star product). In defin- 
ing the Wigner function for the non-commutative pho- 
ton fields, we have to worry about the gauge covariance 
properties. Since the non-commutative photon fields are 
self-interacting much like the usual gluon fields in QCD, 
we define an analogous gauge covariant Wigner function 
as 



d 4 y 



^• fe G«(z)*G^-)(x) (11) 



W ^ k )-J (2 ^ )4 
where we have defined 

G ffi( x ) = U(x,x±)-kF flv (x±)*U(x±,x) 



(12) 



II. WIGNER FUNCTION APPROACH IN 
NON-COMMUTATIVE QED 

Non-commutative QED is described by the Lagrangian 
density of the form 



L=--F ilv *F» v 



(4) 



where 



Ffxv — djj,A v dyA^ j~^[Afj,, A V ] MB (5) 



and the C/'s represent the link operators in the non- 
commutative theory defined along a straight path, 
namely 



U(x,x±) 



(13) 



Here "P" stands for path ordering from left to right and 
it is straightforward to check that under a gauge trans- 
formation, |JSJ, the link operators in (|13fl transform co- 
variantly as 

U (x, x± ) -» n- 1 (x)*u(x,x±)*n(x±) (w) 
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much like in the usual Yang-Mills theory. As a result, 
the Wigner function, (|llfl . transforms covariantly in the 
adjoint representation under a gauge transformation, 

w^{ x ,k) ^n- 1 {x)*w^(x,k)*n{x) (15) 

We note here that, for a link operator defined along a 
straight path, we can also write 

G%\x) = (ct^ D F pv {x)) (16) 

where the covariant derivative is defined to be in the 
adjoint representation. Furthermore, from the definition 
of the Wigner function in l|ll|) , it is easy to check that 

W*„(x,k) = W VM (x,k) (17) 

In non-commutative QED, charge conjugation is defined 
as the simultaneous transformation [22] 

A^x) -► -A„{x), 9^ -9^ (18) 

I 



Under charge conjugation, it is seen, using various iden- 
tities for star products, that 



W pu {x, k) - W c pv {x, k) = W vil (x, -k) (19) 



Note from the definition in l|ll|) that in the limit 9^0, 
the Wigner function for usual QED is charge conjugation 
invariant, while it transforms non-trivially in the non- 
commutative case. 

The derivation of the transport equation for the 
Wigner function, from its definition in Qlljl. follows in 
a straightforward manner. We discuss the details of such 
a derivation in appendix A. Here, we simply note that 
the crucial relations that play an important role in this 
derivation are given by 



D^U(x,x±) = d^U(x,x ± ) - -(A p (x)*U(x,x ± ) - U(x,x ± )* A p (x±)) 



Di x) U(x ± ,x) 



TV 



dt 



(ef 2 ' D F,t V (x)) J *U(x,x±) 



u(x± , x) - — (A„{x ± ) * U(x±, x) - U(x±, x) * A^ (x) ) 



±^y v U(x±,x)*U dt(^*- D F^(x)) 



(20) 



With the relations in (|20() . it is straightforward to show that the Wigner function for the non-commutative photon 
satisfies 



k-DW^(x,k) = ^-j^k? 

2c dk a 



dt (ey dk D F pa {x)^ *W p „(x,k) + W^(x,k)* (J dt (e,'' ,8i '%(i)) 



k" 



(27rft) 4 
d 4 y 



F pa {x)) + (e* 2 F pa (x) 



{2nh) 



e-^ k G$(x)*[j o di((ej 

■*»•* \u(x, x+) * (£><»+>F„ A )(s+) * U{x+,x) * Gy-\x) 
+G { +\x) * U(x, x_) * (a;_) * U(x_,x) 



■G x y\x) 



(21) 



Here dk represents the derivative with respect to k and 
we have used the compact notation introduced in l|ltj|) to 
write the equation in a simpler form. 

Equation l|21|) represents the full transport equation for 



the photon Wigner function in non-commutative QED. 
We note that by using various identities (which we discuss 
in the appendix A) the second term in l|21(l can be written 
as 



4 



ieh d 



d 4 y 



/ dt t{\(ef- D F>>a(x)),U(x ) x + )*(D p F l , x (x + ))*U(x + ,x)\ *G x J-~\x) 

.In l L \ / J MB 



4c dk a J {2nh) 4 

+U(x, x+) * {D p F^(x + )) * U(x+,x) * [(e-'^- D F<^ , G\}-\x) 
(eJ D FP(x)) , G$(x)] mb * U(x, x-) * (D p F x (x-)) * U(x-,x) 



MB 



d'y 
{2irh) 4 



e* r- "^(a;)),l7(a!,a;_)*(Dp^(a;_))*^(x_,a;)l } 

/ J MB J 

^ k [U{x, x+) * [F^(x + ), F$(x + )] MB * U(x + , x) * G x y\x) 



-G$ (x)*U(x,x^)* [F Aff (!_), F m (*_)] MB * U(x- , x) } 



(22) 



In the hard thermal loop approximation, as we have men- 
tioned earlier, it is assumed that the gradients in the 
system are small compared with k/H. Moreover, for the 
semi-classical picture to hold, one also assumes that the 
ensemble average of the covariant derivative, (D W), may 
be considered as being sufficiently small compared with 
(kW)/H. In this approximation, it may be verified that 
the ensemble average of the term l|22[l. namely, of the 
second term in <|21|l is small compared to that of the 



first term. Under the above conditions, the exponen- 
tials exp + (±|y • D) ~ exp^i^lrC^ ■ D) which appear in 
the first term of l|21|) may also be approximated by 1. 
Therefore, if we are only interested in the leading order 
behavior in the hard thermal loop approximation, we can 
neglect the second term in Ij21ll and the classical trans- 
port equation for the Wigner function takes the simple 
form 



k-DW^ix.k) = — — k p \F pa (x)*W llv (x ) k)+W til/ (x > k)*F pa (x) 



d 4 y 



e-^ h G^{x)*F pa {x)*Gy-\x) 



(23) 



Since non-commutative QED is a self-interacting the- 
ory, very much like the conventional QCD, one has to 
consistently separate the gauge field into a background 
part and a quantum part. Normally, this is done by as- 
suming a mean field decomposition of the form 

A^x) = A^x) + a^x) (24) 

where it is assumed that, in this mean field approxima- 
tion, 

(A^x)) = A^x), (a^x)) = (25) 

In making such a decomposition, it is assumed, as is the 
case in the usual background field method, that under a 
gauge transformation, JHJ, 

A^x) -> 9r l {x) * A p (x) *Q,(x) 

+ — n- 1 (x)*d^ i n(x), (26a) 



This is very important in understanding the gauge trans- 
formation properties of the quantities resulting from the 
transport equation. We note that, under such a decom- 
position, 

F^ v {x) = Ffj, v {x) + D^a v (x) - D^a^x) 

1G 

- y c I^W'^WImb ( 27 ) 

where 

D^a u = - — [A^, a„] MB , (28a) 

Ftu, = dfiA-v - d v \ L - ^ [A^Av] MB (28b) 

It is clear that every term in (|27|l transforms covariantly 
under the gauge transformation i|26[l . 
Let us also define 



a^{x) — > VL \x) * dpix) -k fi(x) 



(26b) 



G„ v {x,k) = {W^(x,k))-W^{x,k) (29) 
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where Wfj, v {x,k) represents the Wigner function associ- 
ated with the background field and it is important to 
recognize that 

(W^(x,k)) ^W^{x,k) (30) 
We note that W^(x, k) satisfies the same equation as 



in (1 2 ■" > I to leading order with all background fields. If wc 
now assume, as is conventionally done in QCD, that all 
correlations factors vanish unless they involve G^) (x) or 
Wfivlx, k), then we can write a transport equation for 
Gfiu{x, k) to leading order as 



k ■ DQ^ v {x, k) 



e d 



2c dk a 

- 2 

It also follows from this, that 



k" 



F p<y (x) -kQ pv (x, k) + Q,mu(x, k) * F p<J (x) 



d 4 y 
(2nh) 4 



(^+\x)-kF p<r (x)*G^-\x)) (G (± »W - ^'(x)))]. (31) 



T( Xi k) = - v ^g^(x,k) 



(32) 



would satisfy the equation 



d 



k-DF{x,k) = — -—k p F pa (x)-kJ r (x,k)+J 7 (x,k)*F pa (x) 



2c dk a 
2 



d 'y e -*v-fcJL 
{2ttHY k 2 



{(G$(x)*F pa (x)*G^(-\x)) - (G&(x)^GW(xj)}] (33) 



This equation is manifestly covariant under the gauge 
transformation i|26|) and can be solved order by order in 
powers of "e" to give the ensemble average of the Wigner 
function. 

Given the ensemble average of the Wigner function, we 
can now define a covariantly conserved current as follows. 
First, we note that the current is a four vector which is 
odd under charge conjugation (|18fl . Recalling the prop- 
erties of the Wigner function under charge conjugation, 
(jT§|l . let us define 

J,W - \ j d'keik^^ig^ix^-g^ix-k)) 

= | / d i k6{k )k IJ ,{J'{x,k) -F(x,-k)) (34) 

where it is assumed that we sum over the helicity states. 
We note that this current is manifestly odd under charge 
conjugation. It transforms covariantly under the gauge 
transformation l|26Jl and furthermore, using Q33JI ■ is easily 
seen to be covariantly conserved, 

D^{x) = d,J» - £ [Ap, J*] MB = (35) 



This can, therefore, be defined as the current associated 
with our system and can be determined to any order in 
powers of "e" once J-(x, k) is determined. This is a func- 
tional of Ap(x) and would generate n-point amplitudes 
through functional derivation. These can then be com- 
pared with the leading order calculations from perturba- 
tion theory in the hard thermal loop approximation. 



III. LEADING ORDER AMPLITUDES 

Given the equation for !F{x, k), (|33|l . we can now de- 
termine it order by order in "e" . In this section, we will 
calculate this quantity explicitly to second order which 
would give us photon amplitudes up to the three point 
function for which we have explicit results from the calcu- 
lations in perturbation theory. For simplicity of notation, 
we shall use in the following natural units c = h = 1. 

To zeroth order, can be calculated from its defi- 
nition in using Eqs.. JHJ an d 



T {a \x,k) = 



d 4 y 

(27T)* 



lV ' k p <(Va(z+) - 9xa^(x+)) * d x a»(x-)) 



(36) 



The thermal ensemble average in (|36fl can be calculated in a standard manner using the field decomposition for 
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the quantum field a p . The only important thing to re- 
member is that we want manifest gauge covariance pre- 
served in the calculation. This suggests, as in the back- 
ground field method, that the proper gauge condition on 
the quantum fields should maintain this invariance and, 
in particular, a transverse gauge has to be generalized to 
the form 

D^{x) = (37) 

This would suggest that the polarization tensors, in such 
a case, need not be transverse to the momentum four 
vector and, in fact, can have a longitudinal component 



depending on the background field. Fortunately, up to 
the leading three point amplitudes, this modification does 
not give rise to any contribution and, for all practical 
purposes of our calculations in this paper, we can take 
the polarization to be transverse to the momentum four 
vector. In higher order calculations or in the sub-leading 
terms, however, one has to include such contributions 
carefully. 

With this observation, we note that Eq. 1|36|) can be 
evaluated in a straightforward manner and the thermal 
contribution has the form 



(2?r) 4 J (2tt) 

4 5(k 2 )n B (\k°\), (38) 



(27T) 

where ns(fco) = (exp (ko/T) — l) -1 . 

We can now substitute J-^ into 1)33(1 to determine T^ x \ However, the calculation is not quite as iterative in the 
Wigner function approach, as it normally is in the other way of doing. We still have to evaluate some terms on the 
right hand side of (|33|l that do not involve T directly. In trying to evaluate such terms, we note that, to the leading 
order, we can set the link operators to unity in these terms. This is easily seen from the fact that if we expand the link 
operators to any order in "e" , they will involve powers of y which can be thought of as -8^ acting on the integral. Each 
power of y, therefore, gives a contribution that is more and more sub-leading (since k is large) and, consequently, if 
we are interested only in the leading contributions, we can approximately set the link operators to unity in the second 
group of terms on the right hand side of (|3HJ - With this, it follows that to lowest order in "e", these terms give a 
contribution of the form 



A fc P f d "y c -te* 2 

dk a J (2tt) 4 fc 2 



■'— k " I 77^44e-^' £ -2((^ l a A (x + )-a A a M ( a ; + ))*/ p(T ( a; )*aV(x_)} (39) 



Here, we have identified the Abelian part of the field strength tensor as 

f pa (as) = d p A a (x) - d a A p (x) (40) 

The thermal ensemble average in (|39H can be evaluated in a straightforward manner and the temperature dependent 
part has the form 

4e d 



(2tt) 3 dk a 



k p (6(k 2 )n B (\k°\)f pa (x + 8k)) (41) 



Here, we have used the star product identity With the determination of the second term on the right 

., .. hand side of (1331) to the lowest order, we can now deter- 

e* k - x * f(x) * e lk - x = f(x + 6k) (42) minc ^(i) from as 

with the identification 

(Qky = e^ku (43) 
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k-d^ 1) (x,k)-ie k-A,T {0) 



d 



2(2tt) 3 dk a 



k p (f pa {x) *^(°) +^(°) - 8<5(fc 2 )n B (|fc |)/ pff (x + (44) 



Since is independent of x, its Moyal bracket with k ■ A{x) vanishes and using the form for JF<°) in J3BJ), we 
determine 



T (1) {x,k) = 



4e 1 9 



(2tt) 3 fc ■ 9 afc a 



(5(fc 2 )n B (|fc°|)F(/ pff (x) - /^(x + 9k)) 



(45) 



Before going onto calculate , let us discuss some of 
the features of the results in l|38H and H45|) . We note that 
!F^(x, k) is manifestly covariant under a gauge trans- 
formation, l|2t)l) . However, even though we start from a 
gauge covariant equation, (J2SJ), we see signs of violation 
of gauge covariance in the calculation of (x, k). This 
is manifest more clearly in the right hand side of (|41|l . 
Namely, we note that 

e* k - x *f pa (x)*e~* k - x (46) 

does not transform under a gauge transformation, 1)2(1 [I. 
as we would expect. In a non-commutative gauge the- 
ory, space-time translations form a subgroup of the gauge 
group and, in particular, because of relations like (|42|l . 
gauge covariance appears to be violated. We want to 
emphasize that gauge covariance is manifest before tak- 
ing the ensemble average. However, the naive mean field 
ensemble average seems to be incompatible with gauge 
covariance. It is worth pointing out that this is not a 
problem in the usual QCD where factors such as e ik ' x 
are ordinary functions. This, therefore, is a very spe- 
cial feature of the non-commutative nature of the theory 
and implies that, in such theories, the naive mean field 
ensemble average must be modified. 

For lack of a more fundamental understanding of the 
mean field average method in such theories, we proceed as 
follows. We note from (|42|l that the ordinary plane wave 
function in non-commutative theories leads to a trans- 
lation which does not commute with gauge transforma- 
tion. The simplest way to covariantize such an expression 
would be to replace the coordinate in the exponent of the 
plane wave by a covariant coordinate which would gener- 
ate a covariant translation. Such a covariant coordinate 
can, in fact, be uniquely determined to the leading order, 



from a few general conditions that we discuss in detail in 
appendix B. For the present, however, we simply note 
that we can uniquely identify the covariant coordinate 
with 

X» = x< J - + e9^A v (x) (47) 

where we identify (we will discuss this point more in de- 
tail in appendix B) 

Mx) = A^x) + -Lj^aOfc" (48) 

Thus, we see that the simplest way to covariantize the 
mean field calculations in non-commutative theories is 
to replace 

e ±ik-x _^ e ±i(fe-x+efcxA(a)) /^gs 

where we have used the standard notation of non- 
commutative theories, 

A x B = O^A^Bv (50) 

We note that such a covariantization vanishes in the 
usual Yang-Mills theories simply because 0^ u — 0, but 
is crucial for the covariantization of the results in a non- 
commutative theory. 

The covariantization in l|49|l contributes only at higher 
orders in "e". Therefore, J-^(x, k) is unaffected by this. 
However, in calculating T^ 2 \x, k), we have to take into 
account contributions coming from this in order to main- 
tain covariance under a gauge transformation. At order 
"e 2 ", the leading order contributions from the second 
term on the right hand side of (|33|l have the forms 



2e—k" [ ^-e-^ — 
le dk„ J (27r)4 e P 

x {{(d p a\(x + ) - dxa^ix+f) * f pa (x) -k d x a tJl (x-)) CO v 

-ie {(d fJ ,a\(x + ) - d x a ll (x + )) * [A p (x), A a (x)] MB *d x a M (x-)) 

-ie (([Ap,a\(x+)] MB - [A x , a M (ss+)] MB ) * f p „(x) * «9V(x_)) 

-ie ((d M ax(x+) - d x a fl (x + )) * f pa (x) * [A x , a^(x_)] MB }} (51) 
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Here, (■ • ■ } cov stands for the linear terms in A coming 
from the covariantization discussed in l|49|l . The other 
terms are already of order "e 2 " so that the covarianti- 



zation does not contribute in such terms at this order. 
Equation (|51|) can be evaluated in a straightforward man- 
ner to give 



Aie 2 d 
(2tt) 3 dk a 



k" [S(k 2 )n B (\k°\) 



k-d 



k- (A{x + 9k) - A(x)),f pa (x + 9k) 



MB 



+ [A p (x + 9 k), A a {x + 9k)} MB \ + 



(52) 

where represent terms that are of sub-leading order. Substituting this into . we can now determine to 

leading order, 



f {2) (x,k) = ie—— k- A(x),F (1 \x,k) 
k ■ a L 



MB 



Aie 2 1 d 
(2tt) 3 k-ddk a 
1 



kp{s(k 2 )n B (\k°\) ( [A p {x),A a {xj\ MB - [A p (x + 9 k), A a {x + 9 k)} 



MB 



k-d 



k-(A{x)-A(x + 9k))J pa (x + 9k) 



MB 



} 



(53) 



r 



This is all we need to determine the leading amplitudes 
up to the three point function in the hard thermal loop 
approximation. However, let us first note some of the 
basic features of these results. We note from Ij45f) and 
(|53|l that these Wigner functions have the dipole struc- 
ture characteristic of non-commutative theories. As we 
have alluded to in the introduction, the constituents of 
non-commutative theories can be thought of as extended 



particles and have a dipole structure in the charged sec- 
tor. This is basically reflected in these calculations of the 
Wigner functions. 

Given J 7 ^ 1 ' (x,k),i = 0, 1, 2, we can now construct the 
current up to third order in the coupling constant from 
the definition in lj34J) ■ In momentum space, they have 
the explicit forms 



4 0) H>i) = o 



(54a) 



4 2) (-pi) 



jjiH-Pi) = -8e 2 y dKn B {k°)(l - cos (pi x k))L^(pi, k)L v °{p u k)A v {-pi) (54b) 

16ze 3 J d A p 2 d A p 3 dK n B (k°)S( Pl + p 2 + p 3 ) sin (j^^j 

x (1 -cos(p 3 x kj) ( L^ipx.k)^- + L va (p 3 ,k)^-] A»{p 2 )L Xa {p 3l k)A x {p 3 ) 
\ Pi ' K pi ■ K / 



+ (1 - cos( Pl x k)) k ' A{jP ^> L ^fa, k)L va (p x , k)A u {p 2 ) 
Pi - k 

+ (cos( Pl x k) -cos(p 3 x t f"' A{ P 2) L p<J ( Pl ,k)L^(p 3l k)A x ( P3 ) 

pi - k P2 ■ k 



(54c) 



where we have defined 



dK 



d 4 k 
(2tt) 3 



9(k°)6(k 2 ) 



(55a) 



Ljj,v\Pi k) — T\ pv 



k p p v 
p ■ k 



(55b) 
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It is interesting to note here that even though the Wigner 
functions, J-, have a dipole structure, the odd nature of 
the current under charge conjugation has led to a mani- 
fest quadrupole structure for the currents. This was al- 
ready noted in our earlier paper and here we see explicitly 



the reason for this. 

The currents are functionals of the gauge fields and, 
therefore, by successively differentiating the current with 
respect to gauge fields, we can determine various ampli- 
tudes. For example, the two point function follows from 



IV (p) = 



where we have defined 



sm-p) 



5A U (p) 



A=0 



5A U (p) 



8e 2 

W) 



dKn B (k°) (1 - cos(p x fc)) G M „(p; fc) 



a inUK /,-, = „„, - Kp : + k J p » + = L^k)L;{ P ,k) 



(k-p) {k-p) 
Similarly, the leading order three point amplitude is obtained to be 



6A»( P2 )5A*(p 3 
16 ie 3 



<5 2 4 2) (-pi) 

A=o 5A»(p 2 )6A*(p 3 ) 



sm 



Pi X p 2 



dK n,B 



(*°) {[1 - cos(pi x k)] k x ( Pl ; k) 
k ■ v-i L 



(2?r) 3 V 2 J J ' ' H-Pi 

+ [1 - cos(p 3 x fc)] [hp G v \(p 3 ; k) + k v G^ipi; k) G a x {p 3 ] fc)] 

+ [cos(pi x fc) - cos(p 3 x fc)] ^— — k v G^ipi; fc) G a x (p 3 ] fc) - (p 2 <-> p 3 \v <-> A) 

fc • J>2 



(56) 



(57) 



(58) 



These amplitudes agree completely with the leading 
order perturbative results, arising from the gauge and 
the ghost loops, in the hard thermal loop approximation 
discussed in \2(\ . They are gauge independent and satisfy 
simple Ward identities. For example, from the structure 
in l|57|) . one can easily verify the transversality of the 



photon self-energy 

P M n /w/ (p)=0 J (59) 

Similarly, the identity relating the two and the three 
point functions 



P 3 Vx(pi,p 2 ,P3) = 2zesin I Pl - P2 ) PV(pi) - IV(p 2 )] • KHl'i 



can also be seen to hold. 



IV. SUMMARY AND DISCUSSIONS 



In this paper, we have studied the Wigner function 
approach to the transport equation for photon in non- 
commutative QED. This is complementary to our earlier 
work j2(J in that here we have tried to give a theoret- 
ical derivation for the transport equation. While non- 
commutative QED has many similarities with the con- 
ventional QCD, some differences arise in this method. 
In particular, we have shown that because the gauge 
symmetry in non-commutative QED is intermixed with 
space-time translations, the naive mean field ensemble 



average breaks down and needs modification. While we 
have proposed a covariantization with an interest in the 
leading order hard thermal loop calculations, a system- 
atic study of this question remains to be carried out. 
We have shown that the leading order amplitudes result- 
ing from the current in the Wigner function transport 
equation approach in non-commutative QED agree com- 
pletely with the explicit perturbation calculations up to 
the three point amplitudes. 

We can also derive the force law for a neutral non- 
commutative particle from our calculations. Let us recall 
that in the usual description of the transport equation 
(using dynamical equations for the particles), the current 
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is defined as 

J^x)=efdKk ft f{x,k) (61) 

where f(x, k) represents the distribution function. Com- 
paring with our definition of the current in l|34|). we see 
that 

f(x,k)~±{?(x,k)-?(x,-k)). (62) 

We may now relate the form of the transport equation 
(I33[) with the structure based on the dynamical equations 
(121) for particles in non-commutative QED. As shown in 
[26j , to leading order, such a transport equation would 
be of the form 

k-Df(x,k)+e^(n B (\k°\)X°) = 0. (63) 

Comparing (|63|l with (|33|l , and using the identification in 
(|62p. as well as the relations (|45|l and (|53|l . we find that 

X' 7 = 2 (l - cos k x (ID)) F^fcp (64) 

This force equation for the charge neutral particle, 
agrees completely with what we had obtained in our ear- 
lier paper. As we had emphasized in that paper, we 
see that the nature of the force on such a constituent 
is that of a quadrupole. However, this theoretical deriva- 
tion clarifies how dipole structures inherent in a non- 
commutative theories give rise to a quadrupole effect. 
To understand this, we remark that up to corrections of 
order e, the above result may be written in the form 

X° ~ {k p [ft>{x) ~ r"(x + 9k)]}- {k -> -k}. (65) 



We see that each of the above force terms indeed has 
a dipole structure. However, the charge conjugation 
odd nature of the current converts the dipoles into a 
quadrupole which effectively has the structure of a pair 
of dipoles back to back. 
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APPENDIX A: DERIVATION OF THE 
TRANSPORT EQUATION 

In this appendix, we will give some details on the 
derivation of the transport equation in (|21|l . First, we 
note that, for a link operator along a straight path, we 
can write, 

U(x,x ± ) = pLpJo 1 ^^^ 

= ^L^Sodt^A{x±{i-t)%)\ ( A1) 

where we have identified 

e(t) = ±(i-t)^ (A2) 

It follows now from its definition that, under a general 
change of the end points, the link operator will change 
as 



5U(x, x ± ) = —J dt U(x, x + * S (£(t) ■ A(x + £(t))J * U{x + £(t), x ± ) 
With some algebraic manipulations involving integration by parts, the change can be rewritten in the form 
5U(x, x±) = Y 5xtl ( A f*i x ) * u ( x > x ±) ~ u i x ' x ±) * a m{ x ±)) T 2M SyfJ,U ^ x > x ±) * A ^ x ±) 



(A3) 



l 4- I dts{x+zy 

nc J 



dt 



dt 



-i v U(x, x + 0* (d„A v (x + 0- duA^x + 0) * U(x + £, x ± ) 
Let us next recall that, by definition, 



(A4) 



U(x,x + Z(t)) = P ej 



ft dt' £(t')-A(x+£(t')) 



U(x + £(t),x±) = P [el 



H/ 'd4'f(t').A(x+?(t')) 



(A5) 
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It follows from this that 



dU(x, x + £(*)) _ ie 



dt 



He 



U{x,x + S{t))*£(t)-A(x + £{t)) 



dU(x + £(t),x±) ie- enw^TTf , tn\ \ 



(A6) 



Using these relations, the variation in 1A4|1 can be written in the simpler form 

6U(x,x±) = —5x> i {A tl ,(x)*U{x,x±) - U(x,x±) * A M (x±)) =F -^-Sy^Uix, x±) *A M (x±) 
he 2nc 



ie 



he J 



dt e(t)8(x + O v U(x, x + £(t)) * F fiv (x + £(t)) * U(x + £(*), x ± ) 



It now follows from this that, for Sy^ = and an infinitesimal translation of x^, we have 



D\; ] n,-.,- i: T^y v [j dt(e^ D F^tx))}*U{x,x±) 



(A7) 



(A8) 



where the covariant derivative of the link operator is defined in (|2U[1 and we have used the group properties of the link 
operators in writing the expression in the final form. The second relation in l|20[) can also be derived in an exactly 
similar manner. Given these relations, the transport equation H21|) follows from the observation that 

D { p ] G { ±){x) = (p^U{x,x ± j) *F^(x ± )*U(x ± ,x)+U(x,x ± )*(D^F flv (x ± ))*U(x ± ,x) 

+U(x, x±) * F^y (x ± ) * (p¥>U(x±,x?) (A9) 

The simplification of the second term in H21J1 presented in (|22|l can be obtained as follows. First, we note that the 
factor k p can be taken inside the integral as 



k» 



_d_ 

dy P 



(A10) 



Furthermore, the y-derivatives can be converted to covariant derivatives and, from i|A7|l . we see that we can write 
DMu(x,x±) = dMu(x,x±)±^-U(x,x±)*A fl (x±) = 

D^U(x ± ,x) = d^U(x ± ,x)T ^rM*±)*U{x±,x) = ^-y v U(x ± ,x)* ^ dtt (ef^ D F^(x)) } (All) 



Using these as well as the identity (following from the 
Bianchi identity as well as the equations of motion), 



D 2 F IMU = 2-[F tlX ,Fi,] 



I MB 



(A12) 



one obtains the expression given in Eq. (|22(l . 



APPENDIX B: SOME USEFUL RELATIONS 

In this appendix, we will derive some other useful re- 
lations. First, let us note from iJSJ that 

or, [-i6^x»,f(x)] MB = (OJ) (Bl) 



where we are assuming that is invertible. Using this, 
we see that we can write a covariant derivative (in the 
adjoint representation) with a gauge connection A^ as 



(Daf) = 



—i 



MB 



e 

he 



MB 



-iO-Hx v + ^-6» x A x ),f 



MB 



Here. 



X»=af + — B^AJx) 
he 



(B2) 



(B3) 



is known as the covariant coordinate, since its Moyal 
bracket leads to the covariant derivative. Using relation 
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(|B1(1 repeatedly, it follows that 

e* a * f{x) * e~ ik ' x = (e- kxd f) (B4) 
Similarly, using relation (|B2|) repeatedly, it follows that 



i(k-x + ^-kx A) n —i(k-x + ^-kxA) 

'- k * 6 f(xj) = (el D f(x)) 



(B5) 



where we have used the conventional definition (in non- 
commutative field theories) 



(B6) 



Let us next present the arguments that determine the 
form of uniquely. First, it must transform like a gauge 
field and since the component of parallel to fc M does 
not enter into (|B5(I . without loss of generality, we can 
identify 



Thus, if we write 



k ■ A = k ■ A 



A^ — A^ + 



it follows from (|B7|) that 

k ■ Y = 



(B7) 
(B8) 
(B9) 



Since both A^ and A^ transform like gauge connections, 
it follows that must transform covariantly in the ad- 
joint representation under a gauge transformation. If we 
now assume that Y^ has no explicit 9 dependence except 
for the ones coming from the star products (in order to 
be compatible with charge conjugation), the form of Y^ 
is determined to be 



Since k 2 = 0, on dimensional grounds, we can write 



H(k,D,F) = h{k,D,F) 



k ■ D 



■ g{k,D,F) w (Bll) 



where h,g are dimensionless functions. We next make 
an important observation following from the structure of 
thermal theories in the hard thermal loop expansion. It 
is well known that in such a case only singularities of the 
type ( fc .Q)„ arise. Thus, if we are interested in making 
contact with thermal field theories in the hard thermal 
loop approximation, the second structure in (|B11() cannot 
be present. Furthermore, if we write 



h(k,D,F) = J2 



1 



n=0 



(k ■ D) n 



C n (k,D,F) 



(B12) 



with Co = constant, it is clear that the higher order 
terms in the series will be more and more sub-leading by 
power counting. If we are only interested in the leading 
order terms, we can restrict ourselves to the n = term 
in which case we can write 



Yn — C 



k-D ^ 



(B13) 



Finally, the constant Co is determined by requiring that 
the current following from using such a Y^ must be inte- 
grable (so that it can be interpreted as coming from an 
effective action), 



<5 2 r 



SJu 



5A U 5A^8A V 



8J V 
SA^ 



(B14) 



This uniquely determines the constant Co = 1. It is 
interesting that Co = is ruled out by the integrability 
of the current. Thus, we determine uniquely that 



Y tl =H{k,D,F)F ilv k v 



(BIO) 



where H is a gauge covariant function depending on 



J 



A. 1 1 — -An H~ — Fnjik 

k-D M 



which is the result used in (|48|l . 



(B15) 



Finally, let us indicate how one obtains terms linear in A^ from the exponents involving covariant coordinates in 
Q51JI. Using l|B5(l and setting c = h= 1, we note that to linear order in the gauge field, we can write 



i(k-x+ekx A) 



* fpa (■£) * ^* 



i(k-x+ekxA) 



lin 



lin 



-ie 



Y, ~, [{ (* • 9) n -\k -A) + (k- d) n - 2 (k ■ A)(k ■ d) + ■ ■ ■ + (k ■ A)(k ■ 9)™- 1 } , fj (B16) 
' n\ L 1 J J mb 
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This expression is best evaluated in momentum space, where the right hand side takes the form 
= -ie [ d 4 p 2 d 4 p 3 S( Pl +p 2 +p 3 ) ((-2i S m(^^)))(k-A(p 2 ))f ptT (k 3 ) 



oo 1 



x 

n 

n=l 

-2e J d i p 2 d i p 3 6( Pl +p 2 +p 3 ) sin (il^El ) ( /,. . U ,» , ) )./;„„ (,,, ) 
x g 1 (iA-pi) n -(-*fc-p 3 ) B 



^[»! {ik ■ Pl ) - {-ik ■ P3 ) 
2 ie J d^psSip! + P2 + P3 ) sin (^p) ^^/^(ps) (e* Pl - e^*™) (B17) 



I 

where the last parenthesis in (|B17(I becomes a difference of two 

. . cosines, as is manifested in the last term of in 1541 1 . 

f P M = -i(pz p A„(jp 3 )-p 3a A p (p 3 )) (B18) 

When charge conjugation is implemented in the current, 
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